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ABSTRACT 

We compare new observations on the kinematic characteristics of the damped Lya systems against 
results from numerical SPH simulations to test the predictions of hierarchical galaxy formation. This 
exercise is particularly motivated by recent numerical results on the cross-section of damped Lya systems. 
Our analysis focuses on the velocity widths Aw of w 50 low-ion absorption profiles from our sample of 
z > 1.5 damped Lya systems. The results indicate that current numerical simulations fail to match 
the damped Lya observations at high confidence levels (> 99.9%). Although we do not believe that 
our results present an insurmountable challenge to the paradigm of hierarchical cosmology, the damped 
Lya observations suggest that current numerical SPH simulations overlook an integral aspect of galaxy 
formation. 

Subject headings: galaxies: formation, galaxies: high-rcdshift, galaxies: kinematics and dynamics, 
quasars: absorption lines 

1. INTRODUCTION 



With the apparent agreement between recent CMB ex- 
periments (de Bernardis et al. 2001; Lee et al. 2001; Pryke 
et al. 2001), high redshift supernovae surveys (Perlmutter 
et al. 1999; Riess et al. 2001), and several other cosmolog- 
ical diagnostics (e.g. cluster surveys, Bahcall et al. 1999), 
there is growing optimism that adiabatic Cold Dark Mat- 
ter (CDM) cosmology accurately describes our universe on 
linear scales. Over the past two decades researchers have 
built a theoretical framework for evolving the initial adi- 
abatic fluctuations to the non-linear regime (e.g. Peebles 
1980). To examine the formation and evolution of galaxies, 
researchers have introduced semi-analytic formalism (e.g. 
Press & Schechter 1974) and pursued comprehensive nu- 
merical simulations (e.g. Davis et al. 1985). These galaxy 
formation scenarios are constrained and tested against the 
tremendous database of galactic observations at low red- 
shift. The models have achieved reasonable success in 
matching the luminosity functions, number counts, corre- 
lation functions, and TuUy-Fisher relation of low redshift 
galaxies. On the other hand, there are still serious incon- 
sistencies between theoretical prediction and observation 
on small scales including the inner density profiles of galax- 
ies (e.g. de Blok, McGaugh, & Rubin 2001) and the an- 
gular momentum of spiral disks (e.g. Navarro & Steinmetz 
2000). Nevertheless, this research presents the community 
with important insight into the processes of galaxy forma- 
tion. 

Unfortunately, there are few observations on galactic 
scales at high redshift and therefore few constraints on 

^ Visiting Astronomer, W.M. Keck Telescope. The Keck Observa- 
tory is a joint facility of the University of California and the Califor- 
nia Institute of Technology. 



models of protogalaxies within hierarchical cosmology. The 
most significant population of high z protogalaxies identi- 
fied in emission are the Lyman break galaxies (Steidel et 
al. 1996). These rapidly star forming protogalaxies - se- 
lected with broad band filters designed to identify a strong 
decrement in flux below the Lyman break - have provided 
significant insight into star formation in the early uni- 
verse. Observations of the Lyman break clustering prop- 
erties (Giavalisco et al. 1998; Adelberger et al. 1998) and 
luminosity function (Steidel et al. 1999) have provided 
measurements on the protogalactic star formation rate at 
z > 3 (Steidel et al. 1998; Madau et al. 1996), insight into 
the sub-structure of protogalactic halos (Weschsler et al. 
1998; Bullock, Weischsler, & Somerville 2001), and hints 
at chemical evolution (Pettini et al. 2001). Another active 
area of high z protogalactic research derives from the spec- 
tral analysis of distant quasars, i.e., quasar absorption line 
studies. The damped Lya systems, Lya absorption line 
features with HI column densities in excess of 10^°'^ cm~^, 
are widely believed to be associated with protogalaxies in 
the early universe (Wolfe et al. 1986). Because these sys- 
tems are identified in absorption (i.e. weighted by optical 
depth) , they should probe protogalaxies with a broad dis- 
tribution of physical characteristics (e.g. mass, luminosity, 
morphology) . In turn, the analysis of high resolution spec- 
troscopy allows one to explore the chemical abundances 
(Lu et al. 1996; Prochaska & Wolfe 1999, 2001), metal- 
licity evolution (Pettini et al. 1994, 1997; Prochaska & 
Wolfe 2000; Prochaska, Gawiser, & Wolfe 2001), star 
formation rate (Wolfe, Prochaska, & Gawiser 2001), and 
kinematic properties (Prochaska & Wolfe 1997, hereafter 
PW97; Prochaska & Wolfe 1998, hereafter PW98; Wolfe 
& Prochaska 2000a, hereafter WPOO; Wolfe & Prochaska 
2000b). These observations help describe the processes of 
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galaxy formation at early epochs and establish the starting 
point from which modern galaxies evolved. 

In PW97, we examined the kinematic characteristics of 
17 damped Lya systems and tested several protogalac- 
tic models against our observations. Of greatest impor- 
tance, we tested the predictions from the leading galaxy 
formation model of CDM which characterized protogalax- 
ies as centrifugally supported disks within dark matter ha- 
los (e.g. Kauffmann 1996; Mo, Mao, & White 1998). 
We found that this scenario could be ruled out at very 
high confidence level because the CDM mass spectrum at 
z ~ 2 predicts far too many galaxies with circular veloc- 
ity Vc < 100 km to reproduce the observed velocity 
widths Av of the damped Lya metal-lines. Inspired by 
our results, Haehnelt, Steinmetz, & Ranch (1998; here- 
after II98) examined high resolution (« 1 kpc) numeri- 
cal simulations of a small sample of dark matter halos to 
investigate protogalactic velocity fields in the CDM cos- 
mology. In contrast to the semi-analytic models which 
assumed a rotating disk morphology, H98 found the kine- 
matics of protogalaxies were best described by a combina- 
tion of rotation, infall, and random motions. Therefore, 
the typical velocity width in these protogalactic halos is 
Av instead of the ~ ^Vc predicted by a thick, 

rotating disk. Doubling the average Av value brought 
the predictions of CDM into reasonable agreement with 
the damped Lya observations provided the cross-section 
of protogalaxies scales like a{vc) oc u^-^ (Haehnelt, Stein- 
metz, & Rauch 2000). This latter assumption did not 
arise from the numerical simulations but followed from a 
simple scaling argument between radial scale-length, Vc, 
and the cross-section. However, a more recent analysis 
of the cross-section of damped Lya systems in a series of 
SPH numerical simulations indicates (j{vc) oc w^'^^ (Gard- 
ner et al. 2001, hereafter GOl). In this Letter we test the 
implications of these numerical simulations against new 
observations on the kinematic characteristics of high red- 
shift damped Lya systems. With the addition of nearly 25 
systems since PW98, the damped Lya systems place very 
strong constraints on the processes of galaxy formation in 
the early universe. 

2. ANALYSIS 

In terms of the kinematic characteristics of the damped 
Lya systems, the velocity width Av of the low-ion profiles^ 
places the most fundamental constraint on protogalactic 
models. This statistic probes both the gravitational mass 
of the protogalaxy as well as the specific morphology of the 
neutral baryonic gas. Following PW97, we define the Av 
value as the interval encompassing 90% of the total optical 
depth. Provided the analysis is restricted to unsaturated 
low-ion transitions, this statistic is a robust measure of 
the velocity fields probed by the damped Lya sightlines. 
In Figure 2a we present the A^; values from the Zabs > 1-5 
damped Lya systems presented in PW97, PW98, WPOO, 
13 new systems taken from Prochaska et al. (2001), and 
two systems kindly provided by Ellison et al. (2000) and 
Dessagues-Zavadsky et al. (2000). The values range from 
At; « 20 — 350 km s~^ with a median of «85 km/s includ- 
ing a significant number of systems with Av > 120 km s^^. 
In a future paper (Prochaska & Wolfe 2002), we discuss 

^ Metal-line transitions prominent in HI regions, e.g., Si II 1808 



the kinematic properties of the new systems at greater 
length and pursue other aspects of these kinematic obser- 
vations. 

Now consider the predictions of the damped Lya Av dis- 
tribution derived from numerical SPH simulations of the 
CDM cosmology. Because of the competing expense of 
spatial resolution and cosmological volume, we must com- 
bine the results of two numerical simulations: (1) the high 
resolution, small volume analysis of H98 which demon- 
strated Av « and (2) the cross-section analysis of 
a cosmologically relevant volume by GOl which showed 
(t{vc) oc with (3 « 1 — 1.5. We can convolve these two re- 
sults with the mass spectrum nDM{vc) of the dark matter 
halos at the redshift of each damped Lya system. In the 
following, we calculate nDM{vc) following the prescription 
of Jenkins et al. (2001) and limit ourselves to the ACDM 
cosmology (ag = 1.0, = 0.7, a« = 0.3, A = 0.7). The 
most important free parameter in the analysis is the limit- 
ing velocity vum below which a DM halo will not give rise 
to a damped Lya system. The analysis by GOl suggests a 
value viim ~ 60 km s~^ is required to match the observed 
number density of damped Lya systems njjLAiz)- Other 
analyses, however, have suggested that halos with Vc as low 
as 35 km s~^ will maintain significant reservoirs of neutral 
gas (Thoul & Weinberg 1996; Navarro & Steinmetz 1997; 
Kepner, Babul, & Spergel 1997). In the following we con- 
sider a range of vum values and stress the importance of 
this single parameter on our conclusions for the damped 
Lya systems. 




Ve (km/s) 

Fig. 1. — Log t(uc) for damped Lya systems at redshift 2 = 2 
(solid), 3 (dashed) and 4 (dash-dot). The dotted line demarcates the 
viim value below which a dark matter halo is presumed incapable 
of giving rise to a damped Lyo system. The curves demonstrate 
the predicted dominance of low mass halos in the ACDM cosmology 
for the damped Ly« systems. One also notes the evolution toward 
more massive systems at lower redshift. 

We build the theoretical Av distributions by randomly 
drawing 500 dark matter halos at the absorption redshift 
of each damped Lya system weighted by optical depth 
t{z^Vc) = n]jLA{z,Vc) cr{z,Vc) and assuming Av — |wc. 
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For a{z,Vc) oc v^, we interpolate the (3 values presented 
by GOl at z = 2, 3, and 4 for the ACDM cosmology. Fig- 
ure 1 plots logr(z,Dc) at z = 2 (solid), 3 (dashed), and 
4 (dash-dot) for Vc = 25 — 1000 km s^^, where each curve 
has been normalized to have a maximum optical depth of 
unity. The curves reinforce the fundamental tenant of hi- 
erarchical cosmology: the mass spectrum at high redshift 
is dominated by low mass halos which merge in time to 
form more massive structures. Furthermore, the curves in 
Figure 1 emphasize the importance of the vum parame- 
ter. The optical depth will be dominated by galaxies with 
Vc ~ viim and therefore all predictions on the kinematic 
characteristics of the damped Lya systems are sensitive to 
this parameter. 

Figure 2 presents the Aw distributions for the ACDM 
cosmology for vum = (b) 30 km/s, (c) 45 km/s, (d) 60 
km/s, (e) 75 km/s, and (f) 90 km/s. For each theoret- 
ical distribution we performed a two-sided Kolmogorov- 
Smirnov test against the observed damped Lya Av distri- 
bution. The probability Pks that the two distributions 
could have been drawn for the same parent population is 
indicated in Figure 2. This statistic is most sensitive to 
differences in the medians of two distributions and is a con- 
servative measure; low Pks values are indicative of qual- 
itatively poor matches. The results presented in Figure 2 
indicate that models with < 90 km are ruled out 
at > 99.9% c.l. This suggests a serious failure in our cur- 
rent understanding of galaxy formation within the CDM 
cosmology. 

3. DISCUSSION 

The damped Lya systems present a distinct challenge to 
models of galaxy formation within adiabatic cosmology. In 
particular, the observations impose two well-defined, com- 
peting constraints on galaxy formation scenarios: (1) the 
mass spectrum must include low enough mass halos to 
match the observed number density of damped Lya sys- 
tems, riDLAiz), while (2) the Av distribution limits the 
number of low mass halos which can contribute to damped 
systems. The results presented in the previous section in- 
dicate that current results from numerical SPH simula- 
tions cannot satisfy both constraints. Although a model 
with viim ~ 100 km s~^ could reproduce the median of the 
observed Av distribution, this scenario would severely un- 
derpredict uola- Therefore, physical processes which pre- 
dict large vum values without modifying the cross-sectional 
dependence of the dark matter halos cannot simultane- 
ously satisfy the two constraints. At present, we do not 
consider this failure to be a fatal blow to the paradigm of 
hierarchical galaxy formation, but it does require that at 
least one of our theoretical conditions is incorrect. 

One possibility is that GOl incorrectly determined the 
cross-sectional dependence of dark matter halos containing 
damped Lya systems. In their analysis, GOl had to ex- 
trapolate the results from their SPH simulations to mass 
scales beneath their numerical resolution. The authors 
argued that higher resolution simulations would be un- 
likely to modify their results, but processes on sub-kpc 
scales (e.g. supernovae feedback, tidal gas stripping) could 
have profound implications for the physical nature of the 
damped Lya systems. If such processes implied signifi- 
cantly larger cross-sections for more massive halos, then 



the predicted Av distribution would better match the ob- 
served distribution. To maintain an accurate prediction 
for njjLA, however, the simulations would require a vum 
value even greater than the « 60 km s^^ implied by GOl. 
Because this value already exceeds estimates due to the ef- 
fects of photoionization (Navarro & Steinmetz 1997; Kep- 
ner. Babul, & Spcrgcl 1997; Gnedin 2000), one would 
need to introduce other physical processes to eliminate 
neutral gas from within low mass halos. Perhaps super- 
novae feedback also plays a significant role here. 



Another possible solution is that the Av 



relation 



suggested by H98 is incorrect, particularly in the low Vc 
regime. If physical processes like supernovae feedback dou- 
bled the typical velocity width, then the simulations would 
match the median of the observed Aw distribution with 
'>^iim ~ 50kms~^. At the same time, however, the the- 
oretical distribution would significantly underpredict the 
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Fig. 2. — Velocity width (Av) distributions for (a) our sample 
of z > 1.5 damped Lya systems compared against the predictions 
from SPH numerical simulations of the ACDM cosmology. For the 
model predictions, we randomly draw 500 An values at z^bs of each 
damped Ly« system according to the prescription outlined in the 
text. This includes the assumption that no halos with Vc < 
give rise to damped Ly« systems and we consider a series of 
values: (b) 30 km/s, (c) 45 km/s, (d) 60 km/s, (e) 75 km/s, and (f) 
90 km/s. For ti^m < 90 km/s the models are ruled out at very high 
confidence level by the Kolmogorov-Smirnov test. 
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number of damped Lya systems with Av > 120 km s"""^ 
because the optical depth is so heavily weighted to low Vc 
halos. Furthermore, a model driven by non-gravitational 
processes would need to be finely tuned to reproduce the 
line-profile asymmetries and the correlations observed for 
the low-ion and C IV gas in damped Lya systems (WPOO) . 
These issues could be accurately addressed in a relatively 
inexpensive simulation similar to H98 but including the 
effects of supernovae feedback. We stress that if feedback 
processes are central to an understanding of the damped 
Lya systems then they are likely to dramatically impact 
the formation and evolution of all galaxies. 

The only remaining ingredient of our treatment is the 
mass spectrum of dark matter halos. In this Letter, we 
have adopted the Jenkins et al. (2001) formalism which is 
an analytic fit to the results from the Hubble Volume sim- 
ulations (Evrard et al. 2001). This formalism - at least 
in part avoids the overestimate of low mass halos iden- 
tified in the Press-Schecter formalism. We warn, however, 
that our analysis considers dark matter halos below the 
mass resolution of the simulations analysed by Jenkins et 
al. (2001). It is possible that even this analytic formalism 
overpredicts the mass density of these low mass halos but 
we consider it unlikely that the difference would be large 
enough for the simulations to match the At; distribution 
with viim < 80 km s~^. 

In contrast to the numerical simulations, semi-analytic 
models of galaxy formation can simultaneously account for 
the constraints imposed by the damped Lya observations. 
As described by Mailer et al. (2001), the Av distribution of 
the damped Lya systems can be reproduced by a scenario 
similar to HQS provided a significant fraction of the satel- 
lites residing in dark matter halos have large gaseous ex- 
tent (i.e. cross-sectional dependence a{vc) oc v"^-^). There 
are difficulties associated with this scenario, however. For 
example, these satellites have very low surface density 
which leads to a prediction for the iV(HI) frequency dis- 
tribution, /(-/V), considerably steeper than observed. Fur- 
thermore, the gas distributions adopted by Mailer et al. 
(2001) were driven by the observations and may be diffi- 
cult to motivate physically, particularly if the numerical 
simulations indicate very different gas profiles. Neverthe- 
less, the semi-analytic models may be revealing the path 
that the numerical simulations must ultimately travel; in- 
deed, this is the strength of semi-analytic modeling. 

In closing, we stress the need for a simulation with 
large cosmological volume (e.g. > 20h~^ Mpc) and suffi- 
ciently high spatial resolution to self-consistently examine 
the cross-sectional dependence of the damped Lya systems 
and probe the underlying velocity fields. If the results from 
H98 and GOl stand, the damped Lya observations demon- 
strate that processes critical to the galaxy formation have 
been overlooked by current scenarios. 

We thank E. Gawiser for comments and criticism. J.X.P. 
acknowledges support from a Carnegie postdoctoral fel- 
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